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We calculate the local Green function for a quantum-mechanical particle with hopping between nearest 
and next-nearest neighbors on the Bethe lattice, where the on-site energies may alternate on sublattices. 
For infinite connectivity the renormalized perturbation expansion is carried out by counting all non-self- 
intersecting paths, leading to an implicit equation for the local Green function. By integrating out branches 
of the Bethe lattice the same equation is obtained from a path integral approach for the partition function. 
This also provides the local Green function for finite connectivity. Finally, a recently developed topological 
approach is extended to derive an operator identity which maps the problem onto the case of only nearest- 
neighbor hopping. We find in particular that hopping between next-nearest neighbors leads to an asymmetric 
spectrum with additional van-Hove singularities. 



1 Introduction 

The Bethe lattice plays an important role in statistical mechanics and condensed matter theory. It is defined 
as an infinite tree graph in which each vertex has Z edges, such that any two vertices are connected 
by only one shortest path, as shown in Fig. 1. Several physical problems involving interactions and/or 
disorder can be solved exactly for the Bethe lattice due to its recursive structure, e.g., Ising models [1,2], or 
Anderson localization [3-6]. Furthermore the Bethe lattice is useful as a model for the electronic structure 
of amorphous solids [7]; see Ref. [8] for a recent application. 

In this article we study the spectra of tight-binding Hamiltonians with hopping between nearest (NN) 
and next-nearest neighbors (NNN) on the Bethe lattice. In this context the standard methods of solid-state 
physics which are based on Bloch's theorem cannot be directly applied. This has led to the development 
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Fig. 1 Part of the Bethe lattice with coordination 
number Z = 4. Any two sites are connected by a 
unique shortest path of bonds. Starting from the site 
marked by the open circle, horizontally shaded circles 
can be reached by one lattice step (NN), vertically 
shaded circles by two lattice steps (NNN), and dou- 
bly shaded circles by three lattice steps. Note that the 
lattice is infinite and bipartite. 

of several alternative approaches to calculate the tight-binding spectrum for NN hopping on the Bethe 
lattice [9-14]. However, some of these methods become very complicated for hopping beyond NN due 
to the proliferation of hopping paths. Below we apply and further develop several different calculational 
schemes to determine the Green function in the presence of both NN and NNN hopping. 

In terms of the quantum-mechanical single-particle operator which removes a particle from site 

j and recreates it at site i, the general, translationally invariant hopping Hamiltonian on the Bethe lattice 
has the form 

ij d>0 dij—d 

Here hopping processes between two sites i and j are classified according to their topological distance dij, 
i.e., the number of nearest-neighbor steps of the shortest path between i and j. We will also allow for an 
alternating on-site energy ej, 

i ^ 

where e{ depends only on the sublattice 7 = A, B of the bipartite Bethe lattice to which i belongs. In cor- 
related systems, e.g., for the Hubbard model, on-site energies (2) may be used to detect antiferromagnetic 
symmetry breaking. 

A well-defined limit for infinite coordination number Z results if the hopping amplitudes and Hamilto- 
nians are scaled according to [15] 

td = J^j2 ' H d = Ydfl ' tdHd = l *dHd , (3) 

where K = Z — 1 is the connectivity. In the limit K — ► 00 dynamical mean-field theory (DMFT) [16-21] 
becomes exact. In particular for the Hubbard model the self-energy is then local in space and may be 
obtained from a single-impurity problem with self-consistency condition [20]. In recent years DMFT 
for the Hubbard model on the Bethe lattice has greatly helped to understand the Mott transition from a 
paramagnetic metal to a paramagnetic insulator at half-filling [16-26]. For the paramagnetic phase to 
be stable the antiferromagnetic low-temperature phase of the Hubbard model needs to be suppressed by 
frustration; only then the Mott transition is observable. This can be achieved by including random hopping 
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beyond NN [20,27-30]. In this case the density of states (DOS) remains semi-elliptic, implying that the 
Mott transition in the paramagnetic phase is unchanged. On the other hand, for nonrandom NNN hopping 
the DOS is usually asymmetric, and this is also the case for the Bethe lattice [14]. Moreover, an asymmetric 
DOS is known to stabilize ferromagnetism away from half-filling [31-33]. 

Here, as a prerequisite for any investigations of frustrated interacting or disordered systems, we will 
study the noninteracting Hamiltonian with (nonrandom) NN and NNN hopping 

H = ffioc + txHt + t 2 H 2 = H loc + tlH x + t* 2 H 2 (4) 

and obtain its local Green function Gi(z), defined for Im z ^ by 

Gij{z) := , Gi{z) := G u (z) , Gf{z) := lim G,(z) , (5) 

paying special attention to the limit K — ► oo. Note that due to translational symmetry of the infinite Bethe 
lattice Gi (z) =: G 1 (z) depends only on the sublattice 7 of i. We recall that for only NN hopping the local 
Green function for sublattice 7 is given by [11] 

1 2K(z — f-) 

g y (z):=(i\ -\i) = —— } -11 h = H loc + t* 1 H 1 , (6) 

z-h (K - l)x+ (A + l)y/x- Atfy/x 

where x — (z — e^)(z — £b) and the square roots are given by their principal branches. 

The derivation of the local Green function for t\-t 2 hopping will proceed as follows. In Sec. 2 we use 
the renormalized perturbation expansion (RPE) [11] to obtain a closed set of equations for G^(z), which 
are also related to the DMFT selfconsistency equations. The RPE method is well-suited for the Bethe 
lattice due to its recursive nature, although the classification of paths for t 2 7^ is rather involved. We 
show how to use the RPE result to establish the asymmetry of the DOS. On the other hand, in Sec. 3 we 
use the many-body path integral approach [34] to derive an exact effective action by a recursive method. 
This also yields closed equations for the local Green functions G 1 (z) for any coordination number Z. 
Furthermore a surprising algebraic relation between Green functions for finite and infinite Z is uncovered. 
Finally, in Sec. 4 a recently developed topological method [14] is extended to include Hi oc . We derive an 
operator identity, valid for any Z, that allows one to express G 7 (z) in terms of the known solutions (6) for 
only NN hopping. Our results for the local Green function are discussed in Sec. 5. A conclusion in Sec. 6 
closes the presentation. 



2 Renormalized perturbation expansion 

2.1 Dressed expansion for the local Green function 

In this section we obtain an equation for the local Green function (5) using the renormalized perturbation 
expansion (RPE) [11,35]. In this approach the Green function Gij(z) for H is obtained in terms of the 
Green function G'° c (z) = %G' oc (z) for H loc as 

X S~ll0C I y^loCi /^»10C I \ ^ y^loC * /^loC 1 /^»loC 1 y^loC /n\ 

ij ~ ~r "i z ij {j j "+" / J ^ i ^ik'^'k z kj<^'j T ' ' ' , (_r,j — , \l ) 

, Z Ci 

k 

where we omit the argument z for the moment. This expansion contains terms in which site indices are 
repeated. In a graphical representation this corresponds to paths in which some lattice sites are "decorated" 
with closed paths [11]. One can omit these decorations at the first site i by replacing G' oc by the full local 
Green function Gi. At the next site k, however, the decorations at site i must not be repeated, leading to the 
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replacement of G)° c by the local Green function with site i removed, i.e., by G% = Gfc| Ci=00 . Repeating 
this process one obtains 



Gi — G' oc + y ' Gi tik g\} tki G' oc + y ' Gi tik g\} tkm G[^ t m i G' oc 



(8) 



km 



where the primed sums are now only over non-self-intersecting paths. The RPE is particularly useful in the 
limit Z — > oo which allows the replacement G\ "' — > G°°. This yields an equation involving only local 
Green functions, 



cr(z)- 1 = z 



E 



t lk Gf{z) t kl + Y^Uk Gf{z) t km G™{z) t ri 



km 



(9) 



For the case of two sublattices with on-site energies e A , £b we thus obtain two coupled equations (7 = A,B; 
A = B, B = A), 



G^{zy 1 = z 



F(G™(z),G°°(z)). 



(10) 



This is a closed system of implicit equations for the local Green functions G^°{z). 

Note that the self-consistency equations of DMFT are essentially contained in Eq. (10). In DMFT the 
self-energy is local [20], E y -(z) = £ 7 (z)<5ij for i £ 7, and the local interacting Green function G^(z) 
is given by the Dyson equation G™'^)" 1 = G 7 °(z — £ 7 (z)) _1 = G 1 (z)~ 1 — £ 7 (z), where Q^{z) is 
the Weiss field of an auxiliary single-site problem for sublattice 7. From Eq. (10) one thus obtains the 
self-consistency equation 



G^z)- 1 = z-e y - F(G™\z), Gf{z)) . 



(ID 



Here spin indices were suppressed for simplicity. Detailed discussions of DMFT self-consistency equations 
can be found in Refs. [20, 14]. 



2.2 RPE for tx-t 2 hopping on the Bethe lattice 

For the remainder of this section we consider the Bethe lattice in the limit Z — > 00. As a standard example, 
we first consider only NN hopping (6). Since no closed loops are possible, the only allowed non-self- 
intersecting paths are visits to one of the Z NN sites which return immediately. This yields the well-known 
equation [11] 

g™(z)- 1 =z-e 1 -t?g™(z), g~(z) = 2 ^1^_ (i 2 = 0) (12) 

where the solution is a special case of Eq. (6). 

We now proceed to the case of ti-t-2 hopping, for which the evaluation of the square bracket in Eq. (9) is 
more involved, since it requires the enumeration [36, 26] of several classes of closed non-self-intersecting 
paths which begin and end at site i. First we note that at each site j (7^ i) on the path it is possible to 
make a detour within the same shell, i.e., to one of the yet unvisited NNN sites k of j with same distance 
di,j = di,k from i, as illustrated in Fig. 2a. In the limit Z — > 00 we can take these detours into account by 
replacing G^ by 



Goo poo 1 /-too 1* /-too I poo 1* poo 4* poo I 7 /1 <?\ 



GJ 
i 2 ^ 7 



and only considering non-self-intersecting paths that change shells in every step. Such shell-changing, non- 
self-intersecting paths are referred to as proper paths from now on. They may be drawn using simplified 
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Fig. 2 (a) Example of intra-shell detours in RPE. A non-self-intersecting path, which starts and ends at site i, may 
contain NNN steps (dotted) to sites with the same distance from i, e.g., from j to k. Open and solid circles indicate 
7 and 7 sublattices. (b) Simplified diagram for the proper path corresponding to (a), (c) Rule (14) for proper paths: 
After visiting k and j, the shaded sites are blocked and one must return towards i. From a NN site of i it may then 
either return immediately to i or leave again with NNN steps, eventually taking a NN step and being forced back to i, 
Dashed arrows indicate forced moves. 



diagrams indicating the visited shells only, see Fig. 2b. Closed proper paths, starting and ending at site i, 
are then governed by the following rule: 

Whenever two neighboring sites j and k have been visited, with j further away (14) 
from i than k, all sites on j's branch further away from i than j cannot be visited. 

This rule, which is illustrated in Fig. 2c, is due to the tree-like structure of the Bethe lattice and the require- 
ments that the paths are non-self-intersecting, involve only NN or NNN steps, and change shells in every 
step. It is not affected by intra-shell detours contained in G™. 

Therefore a sequence of outgoing NNN steps, starting at a site in shell or 1 and terminated by a NN 
step must be followed by a NNN step going inward until shell 1 or 0, respectively, is reached. We are thus 
led to the following classification of proper paths by their outermost shell s max , as shown in Fig. 3 and 
listed in Table 1 : 

(a) An even shell s max is reached by taking NNN steps. After taking one NN step inward the path must 
return towards i until reaching shell 1 due to rule (14). Then the path may turn outward again, finally 
using a NN step (al) inward or (a2) outward to reach i's sublattice and return. The reverse path is also 
possible, but in case (al) may be identical. 

(b) An odd shell s max is reached by taking NNN steps and one NN step outward. After returning to shell 
1 due to rule (14), one may turn outward again in (bl) and (b2) similar to (al) and (a2). The reverse 
path is also possible, but in case (b2) may be identical. 

(c) Whereas (a) and (b) involve at least one NN step, there is also one proper path involving no NN step; 
it visits a NNN site and returns immediately. 

Note that paths that start with a NN step are included in the inverted paths mentioned under (a) and (b). 
The contributions of all proper paths are also collected in Table 1, and their sum is 



00 m — 1 00 m 

f (Lr 7 , l_r 7 ) — [ U y,m,k + Lr 7,m,fcJ + [ U 7,m,k + U ~, 



Mb2) 
T ~f,m,k 

m=l k=0 " " m=0k=Q ' 



Wc) 



tfG^ (l-t|G2°) tfG% 



(15) 
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(al) 



(a2) 



(bl) 



(b2) 



(c) 



Fig. 3 Simplified diagrams for proper paths that contribute to Eq. (15). Only paths which change the shell in each 
step are considered; the effect of intra-shell detours is taken into account by working with dressed Green functions 
G^, G^°. Open and solid circles indicate 7 and 7 sublattices, respectively. Dashed arrows indicate moves that are 
forced due to the rules (14). See Table 1 for details. 



Eqs. (10) and (15) are a set of coupled quartic equations for and G?°. Eq. (15), taken together with 
(11), fully determines the DMFT self-consistency equation for t\-t 2 hopping on the Bethe lattice. This 
result differs [36, 14] from the self-consistency equations employed in Refs. [20,27-30], which apply only 
to random hopping. 

The result (15) can immediately be used to determine the moments M n of the density of states. For the 
case £a = € b = we have the large-z expansion 



de 



y> Mn 



-ImG°°(e + i0). 



(16) 



Multiplying Eq. (10) by Gr°, inserting Eq. (16), and comparing coefficients of powers of z we find 



Mi=0, M 2 =t\ z + t* 2 2 , M 3 = {it\ z +t* 2 2 )t* 2 , M 4 = 2tr + Ylt\ l t 



3tn 4 



(17) 



revealing at once that the DOS p°°(e) is asymmetric for 7^ [26], in contrast to the case of random 
hopping [20,27-30]. 

Clearly the RPE is well-suited for the Bethe lattice because rule (14) represents a strict constraint on 
proper paths. However, already for the case of ti-t 2 hopping the RPE requires some care. Furthermore, 
the enumeration of paths for hopping beyond NNN, or for Z < 00, is likely to be very tedious. Another 
drawback is that the RPE only yields an implicit system of equations for the local Green functions. Below 
we derive explicit expressions for them via a different route. 
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path steps factor reached shell 



(a) Smax = 2m 1. to NNN steps outward {t* 2 ) m {G™) m 2m 

(to>1) 2. 1 NN step inward (^)(G°°) 2m -1 

3. (to - 1) NNN steps inward* (^) m " 1 (G?°) m - 1 1 

4. fc NNN steps outward (k < to - 1) (t*,) fc (G°°) fc 2fc + 1 



(al) 5. 1 NN step outward (*i)(G~) 2fc + 2 

6. (fc + 1) NNN steps inward* {t* 2 ) k+1 (G™) k 

reverse path (if different) (2 — <5fc, TO _i) 

= (2 - <5 fc . m -i)(^) 2 (^) 2 " l+2fc (G-) m + fe+1 (G-)™+ fc 



(a2) 5. 1 NN step inward {t*)(G™) 2fc 

6. fc NNN steps inward* (^(G^)*" 1 

reverse path 2 



(b) s max = 2?™+l 1. to NNN steps outward (^) m (G~) m 2to 

(to>0) 2. 1 NN step outward 2m + 1 

3. to NNN steps inward* (t%) m (Gf) m 1 

4. NNN steps outward (fc < to) (t* 2 ) k {G^) k 2k + 1 



(bl) 5. INN step outward (if fc ^ to) (i*)(G°°) 2fc + 2 

6. (fc + 1) NNN steps inward* {t* 2 ) k+1 {G™) k 

reverse path 2 

= 2(1 - 5 k ^){tl)\t^ 2m+2k+ \G^) m + k+1 {G^) m+k+1 



(b2) 5. 1 NN step inward (^)(<%°) 2fc 

6. fc NNN steps inward* (^(G^ 5 )*" 1 

reverse path (if different) (2 — S k . m ) 

G ( *Z, k = (2 - 5 k . m ){tlf{tl) 2m + 2k {G^) m + k {G^r+ k+1 



(c) s max = 2 1. 1 NNN step outward {t* 2 )(G™) 2 

(no NN steps) 2. 1 NN step inward (t 2 ) 



GW = (^) 2 (G~) 



Table 1 Proper paths contributing to Eq. (15) for a Bethe lattice with Z — > oo. Paths (al) and (a2) reach an even 
outermost shell, while paths (bl) and (b2) reach an odd outermost shell; both contain at least one NN step. On the 
other hand, without NN steps only path (c) is possible. Steps forced by rule (14) are marked by an asterisk (*). The 
reversion of each non-symmetric path yields a factor of 2. 
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3 Path-integral approach 

3.1 Green function as path integral 

In this section we use the standard many-body path-integral approach to the local Green function. For 
non-interacting spinless fermions with Hamiltonian H\ oc + H^ n the action is [34] 

S = f ' Vc ] (r)[(a T -M + eJSij + tij ] c, (r) dr , (18) 

Jo ij 

where Ci(r), Ci(r) are Grassmann variables. Fourier transforming to fermionic Matsubara frequencies, 
oj n = (2n + 1)ttT, we obtain S = Y, n J2i 3 [Uj + (e* - i^n - Ci n c jn , i.e., the functional integral 

factorizes with respect to n. We thus consider a fixed Matsubara frequency, set z = \u> n + /i, and omit the 
index n. The local Green function is then given by 

G l (z) = Jl[V[c^c J ]e- § c l c l / JY[V[c 1iCi ] e~ § , (19) 

i J 

with S = J2i Sioc(«) + | Shop(i,j), S loc (i) =-(z- e^CiCi, S hop (i, j) = Uj {ciC 3 +CjCi). Note that for 
this non-interacting system the many-body Matsubara Green function is independent of the temperature T 
and coincides with the Green function defined in Eq. (5). 



3.2 Decomposition of the action for t±-t 2 hopping on the Bethe lattice 

We now consider t\-t2 hopping on the Bethe lattice as in Eq. (4). In this case the following decomposition 
of the Bethe lattice turns out to be useful. For two NN sites i and j we let Z(i\j) denote the set of NN sites 
of i, but with site j omitted. Furthermore B(i\j) shall denote the sites on all the branches which begin at 
the sites in Z(i\j) and lead away from i. These definitions are illustrated in Fig. 4a. 




Fig. 4 (a) The NN sites of i, except j, are denoted by Z(i\j), marked by horizontal shading. B(i\j) denotes the sites 
on the branches starting on any k £ Z(i\j) and leading away from i, marked by horizontal or double shading. S(i|j) 
[Eq. (20)] involves a trace over shaded sites, (b) Here the shaded sites are involved in the trace appearing in H(j|i) 
[Eq. (20)]. When combined with the sites traced over in a(i\j), i.e., the shaded sites in (a), one obtains a trace over all 
sites except i and j, leading to the expression for the local Green function d in Eq. (23). (c) The shaded sites in (a) 
can also be enumerated by combining B(k\i) for all k £ Z(i\j), leading to Eq. (25) which expresses H(i|j) in terms 
of a product H(fc|i)E(fc'|i) ■ • • . 
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Now consider a partial trace where we integrate out the Grassmann variables for the sites in B(i\j), i.e., 
all the sites connected to i except for those on the branch starting at j. We thus define, for NN sites i and j, 

mj)= f II V i^c k ] e-W, (20) 
keB(i\j) 

s ( i \j)= ^2 l Sl <*( k ) + Shop(k,i) + Shap(k,j)] + i ^2 Shop(k,m), (21) 

k£B(i\j) k,meB(i\j) 

i.e., we include only the part S(i\j) of the action in the trace that involves Grassmann variables on the 
branches B(i\j), As a consequence, the only Grassmann variables remaining in are 5j, c, and c 3 , Cj. 

Because they are connected to other in S'(ilj) only linearly, S(z| j) is of Gaussian form, 

Bm«°*(-(*?h)X-r(«j), ^=(5 (22) 

for i e 7, j £ 7. Here the coefficients in the exponent depend only on the sublattices to which i and j 
belong, due to the translational invariance of the infinite Bethe lattice. These coefficients will be determined 
below. 

The partial trace in Eq. (20) serves two purposes. On the one hand the local Green function (19) can 
now be written as a trace over i and a neighboring site j, together with their respective partial traces and 
the remaining parts of the action, 

/Pfr.Cj.Cj.c,-] E(i\j) e-^W-^b)-Shop(M) 
l{Z) JVfaa^Cj] E(i\j) E(j\i) e-s.oc«-s.ocm-s hU pfe) ' <• i} 

This equals (19) because, when considering the two neighboring sites i and j, E(i\j) contains the contri- 
butions of all the sites "behind" i, whereas E(j\i) contains the contributions of all the sites "behind" j, 
as illustrated in Figs. 4a and 4b, and the remaining action involves only i and j because the hopping is 
only between NN or NNN sites. We note in passing that a generalization to hopping beyond NNN, while 
possible in principle, will become increasingly complicated. Inserting Eq. (22) into (23) gives a Gaussian 
integral which can be performed by completing the square, yielding 

G ( z ) = z^-u^- Vj (24) 

7 (zj — Uj — Vy)(z*/ — Uy — Uy) — (ii + Wy + w y )(ti + w 7 + Wj) 

for the local Green function on sublattice 7, with z 1 = z — e 7 . 

On the other hand one can obtain a functional equation for E(i\j) by decomposing the Bethe lattice in 
different ways. Recall that E(i\j) contains the branches which start at the NN sites k of i, where fc 7^ j. 
We can also move one site away from i and consider E(k\i), which contains the branches which start at fc's 
NN neighbors, except i, as illustrated in Fig. 4c. We recover E(i\j) when we multiply the S(fc|i)'s, include 
pieces of the action which connect the B(k\i) among each other and with i and j, and trace over the sites 
fc. Thus we arrive at 



zm= f n v fo> c ^ n e ^ 

kez(i\j) 

Sloc(k) + S b0 p(k,i) + Shap(k, j) + - ^2 5 hop (fc,fc') 




where i and j are any two NN sites. 
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Now can be determined from the above relation, using its Gaussian form (22). Due to the 

translational symmetry of the infinite Bethe lattice it suffices to consider, say, i = € 7 and j = 1 S 7. 
From Eqs. (22) and (25) we then obtain 



f Z I Z - 

J k=2 \ i=2 



S(0|1) cx j Y[ T^[ck, Cfe] cxp ( —if c v^c - ^(£c; + c 4 £) ~ ^ CjMyCj J , (26) 



where if = Z — 1, £ = (ti + w^)cq + t^c\, £ = (ti + w 7 )co + t 2 c\, and the product and sums are 
over nearest neighbors of site other than site 1. The if x if matrix M in the Gaussian integral and its 
inverse are of the form M,j = 5ij a + (1 — 5y) b and (JVf _1 )jj = 5y c + (1 — <5y) d, respectively, with 
c = (a + (if - 2)b)/q, d = -b/q, and q = a 2 + (if - 2)ab - (if - l)b 2 . Here a = Uj - z 7 and b = t 2 . 
The Gaussian integration in Eq. (26) is performed by completing the square, yielding 

H(0|1) oc cxp [ -Kcov^co - ff ( M ~ 1 )y ] = CX P ( -Kcqv^co - ^£ ) . (27) 

Comparing with Eq. (22) we thus obtain the following system of equations: 

u 7 = if i; 7 + if (ti + TtJ 7 ) (t\ + Wy) G 7 , (28) 

w 1 = Kt 2 (h + W 7 ) G 7 , (29) 

W 7 = if t 2 {t\ + w 7 )G 7 , (30) 

u 7 =ifi|<5 7 , (31) 

with the abbreviation G 7 = [z 7 — it 7 — (if — l)^] -1 - 

3.3 Local Green function for arbitrary connectivity 

For easier discussion of the limit Z — ► 00 we use the scaled hopping parameters of Eq. (3). Putting 

z = z + t\jK and G 7 = G 7 / (1 + t%G-i) = — e 7 ~~ u t) 1 we A 11 ^ a f ter some rearrangement, 



_ ^l 2 ^7 (j; ~ ^2 ^7) ^2 ^7 _ - _ 

(l — ^2 ^7 — ^2 ^7)^ ^ — ^2 



G; 1 , (32) 



W-y = 10 7 = — = — ! — , V 1 = h- . (33) 



if 1 — G 7 — G 7 ^f 1 — G 7 

For given z, ca, £b, Eq. (32) is a system of two symmetric equations for Ga, Gb and we denote the appro- 
priate solution by G 7 = f(z , e 7 , e 7 ); note that this function does not depend explicitly on if. Inserting 
Eqs. (32)-(33) into (24) we can then express the local Green function as 

G 1 {z)- l =G- l -^R{G 1 ,G^), G 7 = /(z + |,e 7 ,e 7 ), (34) 

i?(G 7 ,G 7 ) = ^j 2 G 7 (l-t^G 7 ) _^ + (35) 

(1 — ^2 ^~*7)^(1 — P^2 ^7) 1 — ^2 ^7 



1 1 



where again p = Z/K. Taking the limit K — > oo in Eq. (34) we find that G^°(z) = f(z, e 7 , e 7 ), and 

also lmiR-^oo G-y(z) = G^°(z) [Eq. (13)]. This leads us to the remarkable conclusion that the local Green 
function G~ l {z) for arbitrary K is a rational function of the local Green functions Gy(z) for infinite K, 

-l 

(36) 

The function is determined by Eq. (32). However, this is clearly the same implicit equation that 
was obtained in (10) and (15), thus confirming the RPE calculation of the previous section. The function 
G^° (z) is obtained in the next section. 

Finally we note that the unexpected relation (36) can be checked for only NN hopping, 

-l 

(ta = 0) (37) 
which is indeed fulfilled by the results (6) for this case. 



G 7 (z) 



* ; ),G°°(z ' 



K 



K 



>) 



4 Topological approach 

4.1 Operator identities for hopping Hamiltonians 

Recently a topological approach to the tight-binding spectrum of H^ m on the Bethe lattice was developed 
[14], which we will extend to the case of additional A-B on-site energies Hi oc here. 

We begin with a short review of the method of Ref. [14]. Unlike crystal lattices, the (infinite) Bethe 
lattice has the property that the number of paths between two lattice sites i and j which consist of n NN 
steps depends only on the topological distance dij, but not on the relative orientation of i and j. This 
"distance regularity" entails polynomial relations [37] among the tight-binding Hamiltonians Hd- For the 
Bethe lattice they are given by [14] 

L«/2J fj 

E/ n (#i/2) = . H d = U d (H 1 /2)- — U d - 2 (Hi/2), {d > 2) (38) 

where U n (x) are the Chebychev polynomials of the second kind [38]. Similar relations involving the 
Hermite polynomials hold for the infinite-dimensional hypercubic lattice [26,39]. By contrast, Eq. (38) 
is valid for any K, including the one-dimensional chain (K = 1) as well as the limit K — > oo. These 
relations can also be expressed by means of the generating function 

It- 2 °° 1 - t 2 I K °° 

l-xH 1+ Kx 2 ^ l-xH l+ x* ^ 

From these operator identities one concludes that for the Bethe lattice the eigenstates of any hopping 
Hamiltonian (1) are the same as those of H\, and its eigenvalues e(A) can be expressed as a function of 
the eigenvalues A of H\, The calculation of spectral properties, such as the density of states [14] from this 
effective dispersion e(A) is then straightforward. The method works well for arbitrary hopping t d since no 
explicit enumerations are required. 

We now incorporate the effect of alternating on-site energies [Eq. (2)], 

#ioc= £ ^m\ = e + eV, V = ^(-iy\i){i\, fi±^= : jj, (40) 

7=A,B i£7 i 
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where (— 1)* = ±1 for i € A,B, e.g., (— l) 1 = (— l) di -i for some fixed site j £ A. Since the operators 
Hid (H2d+i) connect the same (different) sublattices they commute (anticommute) with V, respectively, 



VH d -(-l) d H d V = 0, 



VH?-(-l) n H?V = 0, 



(41) 



where the second equation follows from the polynomial relations (38) between H d and if™. Together with 
V 2 = 1 we immediately obtain the useful new operator identity 

{aV + PH 1 ) 2n = (a 2 +(3 2 Hf) n , n = 0,1,2,.... (42) 

for arbitrary constants a, (3. This identity makes it possible to reduce resolvent operators involving V and 
Hi to simpler expressions. 

As a simple application let us consider only NN hopping (6). We find by straightforward series expan- 
sion and partial fraction decomposition 



1 



z- e + eV + t\H x 
z-h ~ {z-e) 2 -(e 2 +t* 2 H 2 ] 



It. 



s=±l 



i 



eV 



1 



Vx J s-s/x-tlHi ' 



(43) 



again with x = (z — ep L ){z — e%). We therefore conclude that for NN hopping, surprisingly, any Green 
function for e& ^ can be written as a linear combination of two Green functions for the homogeneous 
case, 6a = 6b, at other arguments. 



4.2 Results for ii-t 2 hopping 

Forii-<2 hopping [Eq. (4)] there are now two possible routes to the local Green function G 1 (z) for arbitrary 
K. On one hand, according to Eq. (36) of the previous section, G 1 (z) is determined by G^(z) alone. The 
latter was already obtained in Ref. [ 14] from the operator identity (39) for the homogeneous case, together 
with a diagonalization of the 2x2 sublattice problem. This gave 

2 



1 



1 



2^ 2(A^ — \\ )tt 



E 



(-1)% - (Af - \)t* 2 ]JX—2^+2 

A,; 



(44a) 



Ai 2 = y]A± \/A 2 - B , A = 1 



(z A - Z B )t* 2 + t\ 

2tf 



B = 



za 



f + 1 



(44b) 



where z 7 = z — e 7 and all square roots are given by their principal branches. While G 1 {z) is given 
explicitly by Eqs. (36) and (44), we note that this approach would be less promising for hopping beyond 
NNN. 

On the other hand we may directly use the operator identities (38) and (42), which provide the relation 

H2 = H\ — p = (vV + Hi) 2 — v 2 — p, with v arbitrary and p = Z/K. This gives us 

t* 2 H 2 =e-t* 2 (p + v 2 ) + tl (vV + Hi) + t* (vV + Hi) 2 , 



H — H\ oc + t x Hi 

with v = e/t\. Performing the partial fraction decomposition for the resolvent we arrive at 

+ 6^ (-1)" 



H £1 



2 

-Y- 



h' 



G^z) 



6 



2 

-Y 



(-!)%(& + £)■ 



(45) 



(46a) 



6,2 = 



-tf ± yftf + Uft* 2 (z - e) + 4if (e 2 + ptf ) 
2t* 



(46b) 



where <? 7 (z) again denotes the local Green function (6) for only NN hopping. We note that this remarkably 
short route to G 1 (z) is straightforward to carry out also for hopping beyond NNN, as well as for off- 
diagonal Green functions. 
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1 ^ 1 ^ 1 — 1 1 1 ^— 1 1 ^ 1 

-2 -1 to 1 2 3-2 -1 oo 1 2 3 

Fig. 5 Local Green function G 7 (w + iO) for the Bethe lattice with ea = eb = 0, t* =0.5, and t| A* = 0, ^ , 1. Left 
column: Z = 4. Right column: Z = oo. Vertical lines mark divergences. 



5 Results for the local Green function 

Using the results of the previous sections, i.e., Eqs. (36) and (44), or Eqs. (46) and (6), the local Green 
function is now available for arbitrary ti-t 2 hopping and on-site energies e A B on the Bethe lattice for finite 
or infinite coordination number. 

In Figs. 5-7 we consider both Z = 4 and Z = oo and compare the unfrustrated case (i| = 0) to weak 
frustration (t^ — and strong frustration (t^ = t*). The homogeneous case, £a = £b, is shown in 

Fig. 5, whereas we chose e A — £b = 2t* in Figs. 6 and 7. 

From these spectra several effects of frustration may be observed. Beginning with the homogeneous 
case (Fig. 5) the imaginary part of the Green function, i.e., the density of states, is no longer symmetric 
if t* 2 7^ 0, as noted in Sec. 2. This is the expected generic behavior for a bipartite lattice with hopping 
between the same sublattices. Furthermore, as discussed already in Ref. [14], a square-root singularity at 
one band edge develops for large enough \t* 2 \- For strong frustration one notices the appearance of several 
additional cusps in both real and imaginary part of the Green function, as well as an increase in bandwidth. 
It is also apparent that in the limit Z = oo the Green function loses some of its features. In this case its 
real part is linear or even flat in part of the band. 

These characteristics persist for the case e A ^ £b (Figs. 6 and 7). In addition the symmetry G 1 (z) = 
Gy(—z) is absent for t\ ^ 0. For strong frustration more van-Hove singularities appear. We also note that 
the band gap, which is present for only NN hopping due to the alternating on-site energies, is closed for 
large enough NNN hopping. 
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Fig. 6 Local Green function G 7 (u>+i0) for the Bethe lattice with Z = 4, ea = — eb = t* = 0.5 and t^/t* = 0, A, 1. 
Left column: Local Green function for sublattice A, Ga(w + iO), which for i| = also appears in Ref. [11]. Right 
column: Local Green function for sublattice B, Gb(w + iO); note the small cusp at u — —0.167. Vertical lines mark 
divergences. 



In summary, the Green function for ti-t 2 hopping shows several qualitatively new properties, which are 
likely to have an impact also on its behavior in interacting many-body or disordered systems, in particular 
for site-diagonal disorder. 



6 Conclusion 

Due to the special topology of the Bethe lattice, the calculation of the Green function of a quantum- 
mechanical particle for hopping beyond nearest neighbors seemed untractable so far. In this paper we pre- 
sented the derivation of an explicit expression for the local Green function for ti~t 2 hopping and sublattice- 
dependent on-site energies e A , e B for arbitrary coordination number Z, employing a set of complementary 
analytical techniques. Implicit equations for were derived by RPE. They also follow from a path 
integral approach, which furthermore yielded the local Green function G 7 for arbitrary Z as a rational 
function of It should be noted that such a functional relation between Green functions at different Z 
is quite unexpected and to our knowledge does not occur for any other lattice. From a topological approach 
explicit expressions for G 1 and were obtained in terms of the Green function g 1 for only NN hopping. 
We found that NNN hopping makes the density of states asymmetric and may induce additional van-Hove 
singularities, increase the total bandwidth, and close gaps that were opened by alternating on-site energies. 
From the experience of these results we conclude that it will be worthwhile to investigate the effects which 
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hopping beyond nearest neighbors has on the physics of interacting many-body and disordered systems on 
the Bethe lattice. 
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